Theory for charge and orbital density-wave states in manganite Lao.sSri 5 Mn04 
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We propose that the high temperature disordered phase in manganites may be described by spin- 
less e 9 electrons. The interactions of spinless electrons may lead to charge and orbital ordered states. 
Our theory explains the Fermi surface probed in angle resolved photoemission spectroscope and the 
observed simultaneous phase transition to charge and orbital ordered state in Lao.sSri.sMnCU. The 
orbital orderings are shown to be d x i_ y i ± d 3z 2_ r 2. 
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Introduction. The manganese oxides are prototype ma- 
terials for the rich physics of the interplay among spin, 
charge and orbital degrees of freedom, which has been 
an important issue in correlated electron systems [U-Q- 
Because of the strong electron-electron interactions, stan- 
dard perturbation theories are not appropriate to study 
these systems. It has been a challenge to develop proper 
theoretical methods to study the phenomena. 

In this paper, we consider single-layered perovskitc 
Lao.5Sri.5Mn04, which shows orbital and charge order- 
ings, followed by spin long range order, as temperature 
T decreases. We propose that the basic physics of the 
high temperature phase and its phase transition may be 
understood in the large Hund's coupling limit, where the 
electronic structure is described by two-fold Mn-3d e g - 
orbital electrons, whose spins are confined to be parallel 
to the local t^g electrons. At the high temperature, the 
local spins are random, and the electronic structure is 
described by the two e g spinless fermions. The transi- 
tion to the charge and orbital ordered states is driven 
by the interactions between these spinless fermions, and 
can be examined by using mean field approximations. 
Our theory predicts a Fermi surface (FS) with a vol- 
ume as twice as large of that of the spin- 1/2 fermions, 
consistent with the recent angle resolved photoemis- 
sion spectroscope (ARPES) and explains the simulta- 
neous orbital and charge orderings in the single layered 
Lao,5Sri,5Mn04. The theory is also consistent with the 
experiments of the bi-layered La1.2Sr1.gMn1.2O7 Q and 
La1.2Sr1.gMn1.2O7 0. 

Manganese Lao.sSri.5Mn04 exhibits two phase tran- 
sitions. The high temperature phase transition occurs 
at T = T co 1=3 220 K, at which the system becomes 
charge and orbital ordered Q. The charge density 
has a checkerboard distribution, and the orbital has an 
ordered wave vector (7r/2,7r/2). The low temperature 
phase transition occurs at T = Tn ~ 110 K, at which 



the antiferromagnetic spin ordering emerges [13 ■ The 
low temperature phase transition has been studied pre- 
viously [llj , here we shall focus on the high temperature 
phase and the related phase transition. 

In Lao.5Sri.5Mn04, the manganese atoms are arranged 
into a two-dimensional square lattice. Each manganese 
atom is surrounded by six oxygen atoms. The electrons 
are well confined in the plane so that this compound is a 
quasi two-dimensional material. The crystal field splits 
the energy levels into a high energy e g orbitals (d x 2_ y 2 
and d 3z 2_ r 2) 1 and a lower energy ti g orbitals (d xy , d yz , 
and d zx ). There are 3.5 d-electrons per Mn atom, which 
suggests 0.5 electron distributed in the two e g orbitals 
whose spin is confined with a local spin- 3/ 2 of the t2 g 
electrons due to the strong Hund's coupling. 

Lai_a;Sri +: ];Mn04 is insulating for all Sr concentra- 
tions x. However, the remnant FS, which is about 190 
meV below the chemical potential, has been probed by 
ARPES [l2j . The observed fermiology consists of a large 
hole-like FS around (7T, 7t) and a very small electron 
pocket around (0,0). The segment of the hole-like FS 
is quite flat, which may induce good FS nesting and lead 
to charge and orbital orderings [12| ■ 

In the large Hund's coupling limit, the spin of the e g 
electron is parallel to the local spins and the charge con- 
figuration on each Mn site i takes the values n, )a = 0, 1, 
with n = c\ a Ci a , and a = 1,2 labeling orbitals d x - y 2 
and d 3z 2_ r 2, respectively. The interaction part of the 
Hamiltonian, including the the Coulomb interactions Uq 
(on-site) and V (nearest neighbor, NN), is given by (see 
supplement material for a detailed derivation), 
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where the sum of i runs over all the sites on a square 



lattice, (ij) represents NN pairs. Note that the spin po- 
larization of the e g electrons is implied. 

The kinetic energy term of the e g electrons can be de- 
scribed by an NN hopping matrix of the two e g -orbitals. 
The single particle part of the Hamiltonian reads, 

H ° = - E cf^/U^+^.c (2) 

<7i is the spin orientation of the t 2g electrons at site i. 
The hopping integrals between the two sites depend on 
the relative spin orientation of the two spins. In the semi- 
classical limit, we have i^fj CT . = cos 13 with 8^ the 
relative angle of the two spins at sites i and j. 

The solution of H$ strongly depends on the spin config- 
urations of the localized t<i g electrons. Here we consider a 
high temperature phase where the spins are random, and 
approximate cos(^y-) s» (cos(f)), which is an averaged 
value of the solid angle and is independent of the pair (ij) . 

Then we have t"^. j a . — (cos (§))£^f, and H is reduced 
to a usual tight-binding model for spinless fermions. The 
pre- factor (cos (f )) represents a reduction of the hopping 
integral due to the random spins. Note that the average 
value of cos (f ) in the solid angle space is 1/3. H then 
can be written as, 

Ho = ~ (cos (°-)) J2 2^7^(fc)4 a %> ( 3 ) 

ka/3 

where t a P is the hopping integral along the x-axis. 711 = 
722 = 7-h 712 = 721 = 7-7 and j±(k) = cosfc^ ± cosk y . 
In what follows, we shall study H = Hq + Hi by solving 
H first and studying the effect of Hj in Eq. (2) from a 
weak coupling approach. 

H can be diagonalized and the eigen-energy is given 

by 

e± = - (cos6/2)(t n +t 22 ) 1+ (k) 

± yV 1 -t 22 ) 2 1+ (k) 2 + 4(t 12 ) 2 7 2. 

The hopping matrix elements are related by Slater- 
Koster formalism [l3l ) if we consider the direct hopping 
between the two NN Mn sites, from which we obtain 
t 22 = i n /3 and t 12 = t 21 = -t n /V3. Hereafter, we will 
take (cos#/2)i n as the energy unit. 

Fig. 1(a) shows the calculated FS for the quarter 
filling of e g electrons, namely 0.5 electron per Mn-site, 
relevant to the single layer Lao.sSri.sMnO^ The area 
inside the FS is one half of the Brillouin zone (BZ), 
consistent with the Luttinger sum rule for the spinless 
fermions. The shape of the FS is in good agreement with 
the ARPES data, which we reproduce here at the upper- 
right of 1/4 of the BZ. Note that the small Fermi pocket 
probed in the ARPES centered at T point may be repro- 
duced in our theory by slightly increasing the number 
of carriers. As we can see, a large segment of the FS is 
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FIG. 1: FS of the spinless fermion model Ho. (a) is for e g elec- 
tron density n=0.5 per site, and (b) is for n=0.6. Red and 
blue colors on the Fermi sheets represent the states mostly 
orbital d x 2_ y 2 or d 3z 2_ r 2, respectively. The upper-right j 
BZ in (a) and (b) show the FS observed in the ARPES ex- 
periments on single layered Lao.sSri.sMnCU [121 and on bi- 
layered La1.2Sr1.sMn1.2O7 [HJ, respectively. The colors used 
in ARPES data represent intensity, while the blue dots in (b) 
are added here to guide the eyes. 

quite flat, and there is a clear nesting at the wave vec- 
tor q = (7r/2,7r/2), which suggests possible instabilities 
toward ordered states. 

Fig. 1(b) shows the FS for electron number 0.6 per 
Mn site, corresponding to the electron density of the bi- 
layer compound La1.2Sr1.sMn1.2O7 Q, where the elec- 
tron layer splitting can be neglected. The shape of the 
Fermi surface is in good agreement with the ARPES re- 
sult of the bi-layered compound, with a Fermi area of 60% 
of the BZ, and the nesting wave vectors cfi — (0.67T, 0) and 
q-2 = (O.671-, O.671-). This is again in support of the spinless 
fermion scenario. 

Orbital density wave instability. We now study the 
effect of Hj. We will first identify the most plausible 
instabilities by using the random phase approximation 
(RPA) analysis. We then apply a mean field approach 
to examine the phase transitions. To study the density- 
wave instabilities, we define the following orbital (o) and 
charge (c) density operators, 

P° = n l+ - = c\ x c l2 + c| 2 Cii, 

Pi = n l+ + m- = c^Cji + c\ 2 c i2 , (4) 

where the orbitals + and — are linear combinations of 
the orbitals d x 2_ y 2 and d 3z 2_,.2, 

4 = -^(4±4). (5) 

As it will become clear later, the orbital ordering in this 
problem is associated with orbitals + and — , instead of 
1 and 2. We introduce a static susceptibility matrix x, 
whose element is defined as 

1 f' 9 

Xaa'^'M) = 2 J dT { T TPaa'{q-,T)p l , li ,(-q^)) . (6) 
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where p aa ,(q) = Y<k c \ + ^ a c k a > 

The orbital and charge susceptibilities are then given 

by 

an, Xaa,ftft 

x°(q) = lEa^Xaa.wCg)- (7) 
Within the RPA, we have 

x = (i + x {0) u c )- 1 x {0) - (8) 

In the above equation, / is an identity operator, and x^ 
is the matrix of the bare susceptibility, 

X2L„(0 = Jf £ C(* + fla£(*K*(*K.(* + fl 

kmn 

x[f(e n (k + q)) - f(e m (k))]/[e m (k) - e n (k + q) + ir,}, (9) 

where m and n are the band indices, and a" n {k) — 
a, k\m, k^ is the orbital weight. We arrange the matrix 
index from 1 to 4 as (a/3) = (11), (22), (12), and (21). 
The interaction matrix U c is of the form U c = U 1 © U 2 , 
where U 1 = V(q)a + (V(q) + U )a u and U 2 = -U ao 
with do an identity matrix and o~\ the first Pauli matrix. 

In the matrix representation described above, the 
upper-left 2 by 2 block in x describes charge part and 
the lower-right block describes the orbital part, as we can 
see from Eqs. ([7]). While U c is block diagonal, x^°Hq) is 
generally not block diagonal, so that the charge and or- 
bital are coupled in the response functions. A special 
case is at q x = ±<7y, where the off-diagonal components 
of x vanishes due to the symmetry in the band struc- 
ture [lj|. This property makes the study of the insta- 
bility at q = (7r/2,7r/2) and q — (tt,tt) simpler. We will 
show that in this case, the inter-orbital nesting connect- 
ing the different segments of the FS between orbital 1 
and orbital 2 favors the ± orbital ordering. To illustrate 
this point, we define orbitals = cos(ip)c\ + sin(t/j)c\ 
and cL = sin( , 0)c{ — cos^)^, with orbital density p°^ — 

c\c a — CpCp. Note that in the vicinity of orbital-density- 
wave instability, the nesting is between electron states 
with mainly orbitals 1 and mainly orbital 2, we have 

X% ~ C0S(V>) Sin(^) Xaa^fc, (10) 
afj, 

which reaches its maximum with i\) — tt /A. Therefore we 
conclude that the ordered orbitals are + and — (for a 
more detailed discussion, see fl6|). 

We have found three types of instabilities in our cal- 
culations, namely the orbital ordering at (tt/2, tt/2) and 
at (tt/2, 0), and the charge ordering at (tt, tt). Note that 
the orbital orderings are related to the FS nesting, while 
the charge ordering is not. In Fig. 2 we plot the suscep- 
tibilities at corresponding wave vectors as functions of 
interaction strengths. As we can see, the orbital suscep- 
tibilities at (tt/2, tt/2) and (7r/2,0) are greatly enhanced 
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FIG. 2: T = RPA susceptibilities, (a): Orbital susceptibil- 
ity X°(<t) at <f = (tt/2, tt/2) and (n/2,0) for orbital ordering 
between + and — . The dotted line is the susceptibility for or- 
bital ordering between orbitals 1 and 2. V — 0. (b): Charge 
susceptibility x c (<l) at q = (tt, 7r) with Uo = 0. 

by the inter-orbital repulsion Uo, and the susceptibility at 
(tt/2, 7t/2) is much larger at large Uq with a critical value 
of Uq ~ 4 for the ordering. Note that the orbital sus- 
ceptibility based on the ordering between orbitals d x i_ y i 
and d3 Z 2_ r 2, X^=o( 7r /^' 7r /2), is much weaker as we can 
see from the dotted line in Fig. 2(a). For the charge 
ordering at (tt, tt), as seen in Fig. 2(b), x c diverges at 
V ~ 1.1, which indicates a phase transition to (tt, tt) 
charge ordering. 

Phase diagram and phase transition. The RPA cal- 
culations above have indicated two possible major in- 
stabilities, the (tt/2, tt/2) orbital order (OO) and (tt, tt) 
charge order (CO). Below we use mean field approach 
to examine the interplay between the two orderings. We 
introduce the two mean fields 

(p°) = (c^ca + c\ 2 cn^j = p D cos(qi ■ r, + 0), 

(Pi) = ( c !i c *i + 4^2) = Pc cos(q 2 • n) + p, (11) 

with qi = (tt/2, tt/2), q 2 = (tt,tt), and p = 0.5. p 
and p c are the order parameters of charge and orbital, 
respectively, while <fi is the phase shift in the real space 
of the orbital order. The mean field Hamiltonian then 
reads 

Hmf = H - Y [Po eI0 ( c k,2 C k+qi,l + c k,i Ck +<u. 2 ) + H.c. 
k 

-A(V-U /S)p c Y,ct a c k+cl2a , (12) 

kn 

The self consistent equations for the mean fields are 

P" e ^ = f £ ( C k+ qi ,l C k,2 + 4+ qi , 2 Ck,l) 
k 

Pc = AT £ (Ck+q 2 ,aCk, Q ) ■ (13) 
ka 

By solving -£/mf together with the self-consistent equa- 
tions (|13[) . we obtain the zero temperature phase dia- 
gram, which is shown in fig. [3] In the calculation, we 
found only two possible phase shift <j) for the orbital 
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FIG. 3: (a). Phase diagram of spinless fermion model at zero 
temperature. CO: charge-ordered phase; OO: orbital-ordered 
phase, (b) Shapes of ordered orbitals + and — . (c) Illustra- 
tion of the ordered states in real space. Electron charge is 
represented by the size of the circle. Orbitals are represented 
by colors: Blue for dominant orbital +, maroon for dominant 
orbital -, and grey for orbital-disordered site. 



ordering, <j) — 7r/4 and <f> = 0, which are denoted as 
00(7r/4) and OO(0), respectively, in the phase diagram. 
The real space modulation of each phase is sketched in 
Fig. Etc). One of the main features of the phase diagram 
is that the system is in the coexistence phase of CO and 
OO(0) in a large parameter space of (Uq, V). The phase 
with only orbital ordered only appears in a tiny phase 
space with very small V and large Uq. We also note that 
there is a sudden change on the orbital ordering phase 
from 00(7r/4) in the absence of CO to OO(0) in the pres- 
ence of CO. Below we shall provide some understanding 
of the later. Let us first consider the orbital ordered only 
phase. The preferred phase 00(7r/4) may be understood 
as the result of losing less kinetic energy due to the orbital 
ordering. The amplitude of the orbital order parameter 
(p°) for the 00(7r/4) phase is p j\[2, while the amplitude 
for the OO(0) phase is p . However, the situation is very 
different in the presence of charge ordering. In that case, 
the local orbital order (p°) is bound by the local charge 
density of electrons (pi). Because (pi) are reduced on 
some sites, the charge ordering suppresses the 00(7r/4) 
phase. On the other hand, the OO(0) is consistent with 
and may even be enhanced by the charge ordering. In the 
limit of strong charge ordering p c = 1/4, the kinetic en- 
ergy term diminishes, and (p a ) — 1 in the phase OO(0), 
in comparison with a maximum value of (p ) = 0.5 in 
the absence of charge ordering. In other words, the pres- 
ence of charge order will induce the 00(7r/4) phase. The 
transition from 00(7r/4) phase to CO+OO(0) phase is 
the first order. 

In Fig. 4(a)-(c), we plot the orbital and charge order 
parameters as functions of V for various Uo at T=0. At 
small Uq — 1, as V increases, CO develops first followed 
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FIG. 4: Panels (a), (b), and (c): V dependence of orbital 
ordering p (blue curves) and charge ordering p c (red curves) 
at U = 1. (a) at U = 0; (b) at U = 3; (c) at Uo = 5.5. 
Panel (d): Temperature dependence of orbital ordering p 
(blue curves) and charge ordering p c (red curves) at Uo = 
4,1/ = 0.5. 

by a coexistent phase with the OO(0) order. At Uo = 3, 
the transition to the charge and orbital ordered state is 
simultaneous as V increases, and is first order with clear 
jumps in the order parameters. At large Uq — 5.5, we 
have only orbital ordering at small V, and co-existent 
phase with charge ordering. And at the charge ordering 
point, the orbital order parameter has a change in both 
the phase (not shown here but discussed before) and its 
magnitude. In Fig. 4(d), we show the order parame- 
ters as functions of temperature for (Uo — 4, V = 1), 
to illustrate the simultaneous first order phase transi- 
tion of the orbital and charge orderings at finite tem- 
perature, which may explain the simultaneous orderings 
observed in experiment of Lao.5Sri.5Mn04 [9(. Note that 
the mean field theory predicts a large first order transi- 
tion, which may change the electronic structure dramati- 
cally under the transition. Since mean field theories usu- 
ally over-estimate the order parameter jumps in the first 
order phase transition, the interpretation of this should 
be more cautious. 

We now discuss the ordered orbital characters of the 
system. Different from the usual rotational invariant 
spin- 1/2 space, the kinetic energy term is not symmetric 
with respect to the rotation in the pseudo-spin e g or- 
bital space. Therefore, there is a selection of specific or- 
bitals for the orbital-density wave ordering. In the early 
theories the ordered orbitals have been predicted to be 
d-3x 2 -z 2 and <^3y 2 -z 2 However, in more recent x-ray 
analysis (l8| . it is found that the ordering of d x 2_ z 2 and 
d y 2_ z 2 is much more compatible with linear dichroism 
measurements. Note that the two orbitals d x 2_ z 2 and 
d y 2_ z 2 are not orthogonal. To further examine this is- 
sue, we have performed the mean field calculations to 
examine the ordering between a general linear combina- 
tion of d 3z 2 and d x a_ y 2, and have found that the order- 
ing between + and — has the lowest energy. Therefore, 
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in contrast to previous arguments that the ordered or- 
bitals are non-orthogonal, we propose that at tempera- 
ture Tjv < T < T co , the ordered orbitals are orthogonal 
orbitals + and — , which are essentially equal mixtures of 
d 3z 2_ r 2 and d x 2_ y i. The shape of each orbital is plotted 
in Fig. 3(b). 

In summary, we have studied the high temperature 
states and the high temperature phase transition in sin- 
gle layered manganite Lao.5Sr1.5Mn.O4. We propose that 
the basic physics of the high temperature states in this 
class of strongly correlated systems may be described by 
an effective spinless Hamiltonian with two e g electrons. 
The FS calculated within the model in the absence of 
the remaining interaction is in good agreement with the 
ARPES data. The orbital and charge density wave or- 
derings may be studied based on the model with the re- 
maining Coulomb interactions. Our mean field theory 
predicts a simultaneous charge and orbital orderings in a 
large parameter space, consistent with the experiments in 
Lao.5Sri.5Mn04. The ordered orbital characters are pre- 
dicted to be a linear combination of d x i_ y i and dj, 7 2_ r 2 
with the equal amplitudes, which may be further justi- 
fied in experiments. The effective Hamiltonian may be 
applied to other systems of strong correlations. 

Supplementary material: effective spinless model for 
manganites with strong Hund's rule coupling 

Wc consider a model Hamiltonian H = Hq + Hi, with 
Hq the single particle part to describe conduction e g 
electrons and Hi the interaction part to describe the 
Coulomb repulsion between two e g electrons and an on- 
site spin coupling of an e g electron and spin-| of the 
localized t 2g electrons. We first consider the interaction 
Hamiltonian, which is given by 

Hi =U ^2 n ia^n iai + (U' - - J) ^2 nnn i2 

ia i 
i i 

ia <i-j> 

(14) 

where U, U' are on-site intra- and inter- orbital direct 
Coulomb repulsive interactions, respectively, and J > 
the exchange Coulomb interaction or the Hund's rule cou- 
pling. By symmetry, U = U' + 2 J. V is the nearest 
neighbor (n.n.) site Coulomb interaction. Si a is the spin 



of an electron of orbital a at site i. We denote a = 1 
for d x 2_ y 2 orbital and a = 2 for d 3r 2_ z 2 orbital. Si is 
the spin- 1 of three localized t2 g electrons at site i, and 
n ia = n ia ^ + n ia ^ is the total electron number opera- 
tor for a given orbital. In our model, the inter-orbital 
Coulomb repulsion between e g and t 2g electrons is a con- 
stant, which can be absorbed into the chemical potential. 

We consider the large Hund's coupling limit, where U, 
U', J are much larger than the kinetic energy term Hq 
below. We shall assume, however, U' — J to be compa- 
rable with the kinetic energy, and may even be treated 
as a perturbation from a technical point of view. We 
may argue for this limit that in a metallic phase, the 
Coulomb interaction U' has a good screening, while the 
Hund's coupling J is not screened, so that U' — J could 
be small. We will further discuss this point in the sum- 
mary of the paper. In this limit, Sj Q is parallel to Si, and 
doubly occupied e g electrons on the same site is allowed 
because it costs an energy of U' — J. Since the local 
spin degrees of freedom of e g electrons are frozen, the e g 
electrons behave like spinless fermions. We remark that 
this scenario is consistent with the recent ARPES mea- 
surements. The average number of e g electrons in the 
single layer compound LaQ.sSri.sMnC^ is 0.5 per Mn- 
site by valence counting. The Fermi volume measured 
in the ARPES is expected to be 25% of the area of the 
first Brillouin zone for spinful fermions and 50% of the 
zone area for spinless fermions according to the Luttinger 
sum rule. The ARPES data for the single- layered man- 
ganites [IH shows that the area enclosed by the Fermi 
surface (Fermi lines in 2D) is about 49% of the zone area, 
in excellent agreement with the spinless fermion scenario, 
but not spinful fermion scenario. Similarly, the ARPES 
data for the bilayer compound La1.2Sr1.gMn1.2O7 give us 
an estimate of the area enclosed by the Fermi surface to 
be about 60% of the Brillouin zone, consistent with the 
spinless fermion scenario for this compound where the 
average number of e g electron is 0.6 per Mn site (3. [T3|. 
Note that the spin degrees of freedom of the e g electron 
is frozen only locally, and the spins at different Mn sites, 
hence the spins of e g electrons at different sites, may have 
different polarizations. Hi in the large Hund's coupling 
limit then takes the form, 

Hi = Up ^ nan i2 + -77 ^ {ni 1 +n l2 )(n jl +n j2 ) 1 (15) 

i <i,j> 

with Uq = V — J, and the spin polarization of the e g 
electrons is implied. 
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